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Abstract 



We use a recently developed action principle in spaces with curvature and 
torsion to derive the Euler equations of motion for a rigid body within the 
body-fixed coordinate system. This serves as an example that the particle 
trajectories in a space with curvature and torsion follow the straightest paths 
(autoparallels), not the shortest paths (geodesies), as commonly believed. 
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1) Introduction 

Since their discovery more than two centuries ago, the Euler equations 

L + flxL = M (1) 

have played a key role in understanding the rotations of a rigid body around a fixed point 
O. (see for example [Q] - and the references therein). The vectors refer to the body-fixed 
system, L being the angular momentum, f2 the instantaneous angular velocity, and M the 
moments of the external forces. The dot denotes differentiation with respect to the time t. 
The unit basis vectors Si(i = 1,2,3) in B may be assumed to point along the eigenvectors 
of the inertia tensor with respect to O. Then the angular moment L has the components 
Li = IjO* (no sum over i). 

To describe also the translational motion of the rigid body one usually chooses the point 
O to coincide with the center of mass moving through space. The motion satisfies the 
additional equations: 

P + fixP = F (2) 

where P = MV is the linear momentum, M the body's total mass, V the velocity of 
the center of mass, and F the driving force. Presently, the equations (^ and (0) play an 
important role in missile dynamics analysis [||]. 

Usually, the two sets of equations are derived as the Newton equations of motion in the 
stationary reference system S: 

1 = m, (3) 
P = f (4) 

Following general conventions, we distinguish vectors and tensors in S and B by using for 
the small letters with Greek indexes and the capital letters with Latin indexes, respectively. 
An transformation 

ft* = e> M , (5) 
V* = e> M , (6) 



carries the of the angular and linear velocities uj^ and v M to the moving body-fixed system 
B (see § - i). 

The 3x3 matrix elements e l M and e\ depend on the coordinates q 11 in the system S. 
They satisfy d[ u e\] ^ 0, d[ u e\] ^ 0, making the transformation (|5|) and ([]) anholonomic 
coordinate transformations (see [[[J- [[7[] for the general theory of the dynamics in anholonomic 
coordinates is presented). 

Under the transformation (f|), the equations (|]) go over into (1) and the additional term 
(17 x L) fc = e k ijVt l U arises as the moment of the gyroscopic force. The additional term 
(f2 x P) fe = e k ijVL l Pi in equation (Q) arise similarly. Both terms are a consequence of the 
anholonomy of the transformation (|jj), (^). 

Within the stationary system S, Hamilton's action principle serves to derive the equations 
of motion (|3]) and (|4]). If one transforms the classical action to the system B, however, 
the description involves nonholonomic coordinates and a naive application of Hamilton's 
principle produces wrong equations of motion lacking the additional gyroscopic moments. 
In 1901, Poincare showed || how to vary an action expressed in terms of nonholonomic 
coordinates. Following his treatment, one certainly recovers the gyroscopic forces. His 
treatment is reviewed in the Textbook 

The purpose of this note is to point out that Poincare's treatment may be viewed geo- 
metrically as an application of a recently proposed action principle in spaces with curvature 
and torsion to the spinning top within the body-fixed reference system. The motivation for 
such a consideration derives from the fact that, in the literature on gravity with curvature 
and torsion || (for the geometry of such spaces see 0), there is a widespread belief that in 
spaces with torsion, spinless particles move on shortest paths. However, it was discovered in 
Ref. |IIJ (when solving the path integral of the hydrogen atom) that the correct trajectories 
are the straightest paths in a given geometry. In Ref. |TI[, a classical action principle was 
found to comply with this physical fact. 

The key to the new action principle is the observation, that a space M q with nonzero 
Riemann curvature and torsion may be mapped locally into a euclidean space M x by an 



anholonomic transformation x l = e % ^(f [fL2 |, [[I0|j . For the particle with mass m in M x , 



the equations of motion x % = yield straight-line trajectories. Under an anholonomic 
transformation, these lines go over into autoparallels satisfying the equation + T a ^q OL q 13 = 
with = ei^d a e L j3 being the Cartan connection. 

The orbits in the euclidean space M x may be derived from Hamilton's principle 5A X = 
applied to the classical action A x = dt^mx l x l . Under the anholonomic transformation 
x l = e\q^, this action goes into the curvilinear form A q = J** 2 dthmg^q^q 11 '. In the space 
M q , a naive application of Hamilton's principle 5A q = produces wrong equations of motion 
in the space M q . One finds q^ + T a ^q OL q 13 = which are the equations for the geodesies 
rather than the autoparallels; they lack torsion force 2mS a ^ $q a (f . 

The problem of describing the dynamics of a rigid body within the body-fixed frame 
may be formulated is somewhat analogous. Here the parameter space of Euler angles has a 
Riemann curvature. By going to the body-fixed frame and describing the system in terms of 
non-holonomic coordinates, the space becomes affine-flat but possesses torsion. The affine 
connection in the space of anholonomic coordinates, is obtained from the spatial derivatives 
of the transformation matrices in (5) and (6). The associated Cartan curvature tensor 
R liU \ K = efid^dy — d u dfj_)e\ vanishes, making the connection affine-flat. The antisymmetric 
part of the connection = — T^ 7 ), which is a tensor, is nonzero giving rise to 

a nonvanishing Riemann curvature tensor R^x 1 ^ 0. The latter is formed from the Levi- 
Cevita connection, also called Christoffel symbol T^ uX = \ (d^g u \ + d u g^\ — <9a5v), where 
5V = e V e ^ i s the Riemann metric in the space of anholonomic coordinates. 



We shall apply the variational principle of Ref. |TTJ and derive, within the body-fixed 
reference system B, both the Euler equations (1) and equations for the translational motion 
(2). 

2) The SO (3) geometry in a rotating body-fixed system produced by an anholo- 
nomic transformation from the rest system 

Consider first a rigid body rotating around a fixed point. Within the body-fixed sys- 
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tern B, we introduce anholonomic coordinates corresponding to the transformation (|5|). 
They are the components of the body's rotation vector in the axis-angle parametrization. 
The anholonomic transformation (^) defines their infinitesimal increments <f$* = Q l dt. 
For a precise specification, let us go to the system S where the standard Euler an- 
gles a = ip 3 , (3 = if 2 , 7 = if 1 [10] parametrize (holonomically) the body's configura- 
tion space M*- 3 ) = SO (3). The components of the angular velocity in the system B are 
Q 1 = Slip 2 — cis 2 ip 3 , fi 2 = Clip 2 + sis 2 ip 3 ,Q 3 = ip 1 + c 2 ip 3 (with the short notation 
c M = cosv? M , s M = simp' 1 ). The basic relation 



d¥ = eV(vW 



(7) 



has coefficients e* M which form the matrix 



si -cis 2 



V 



(8) 



Ci sis 2 

1 c 2 

The symbol d stands for increments do not belong to an integrable function; they do not 
satisfy the Schwarz integrability condition which would read, with the coefficients of (7), 
dy&y - d v e\ = 0. 

The metric g = (g a p) and its inverse g^ 1 = (g a ^) have the matrices 



1 

c 2 1 



' 1 o 



-1 



s 2 2 
-c 2 1 



(9) 



/ 



where \g\ = det g = s 2 2 ; g 1 = \g\ 1 is the determinant of g. From g we calculate the 
Christoffel symbol: 

\g\ 



The associated Riemannian Ricci tensor is R a p = \g a p with a constant nonzero scalar 
Riemannian curvature R = 3/2. The only nonzero components of the affine-flat Cartan 
connection are 



T12 1 — r 2 3 3 — — , r 12 3 — T23 1 — — , r 13 2 — s 2 . (10) 

^2 S2 



From these we find the antisymmetric part is the torsion tensor S a p~f = — |y |g|e a( g 7 . A 
transforming to the local basis e % yields the so-called object of anholonomy M- M 

2S ijk = 2e?e$4S afh (11) 

whose explicit form is 

where denotes the Levi-Civita antisymmetric tensor in Cartesian coordinates (6123 = 1) 
M. This geometry prepares the ground for the action principle to be developed. 
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3) Action principle for the rotational motion in the body-fixed system B 

Pure rotations around an arbitrary fixed point of the body are governed by the Euler equa- 
tions ([[]). The kinetic energy is T = ^IffiQ*. In S, the Lagrangian is Ls = ^g^ip^ip" — U((p) 
with U being the potential energy and = Jje* M e*„ the "kinetic-energy-metric" . The latter 
metric differs from geometric metric g^ v of Eq. (0).[] In this note, the kinetic-energy-metric 
9^u 09 wm pl a y no r °l e - The potential term U(ip) is irrelevant for our considerations and 



will be dropped. 

An application of Hamilton's action principle 8 As = to the classical action As = 
J t * 2 Lsdt in the system S certainly yields the correct equations of motion. Under the trans- 
formation (H), they become the Euler equations (H). 

Let us now transform the action As via (|5]) to the body-fixed system B. The result is 
very simple: 

A B = t dtL B = f 2 dth¥¥. (13) 

By applying naively Hamilton's action principle we would find the equations I$ % = LSl % = 
for each i which are not the correct Euler equations (|l|) — the gyroscopic moments being 
missed. The contradiction is caused by the fact that the variations 8ip^{t) in the space of 
Euler angles M v and the variations in the space of anholonomic coordinates M$, 

are related with each other in a path dependent way ("nonlocal" on the time axis) ||11|| . 
Explicitly, there is the functional equation 

+ f 2 e l Mt)W{t)dt. (14) 

Its most important consequence is that closed paths in the space are not mapped into 
closed paths in the space M$, due to a nonvanishing Burgers vector b l = § e^dip^ ^ 0. 

The usual Hamilton action principle in the system S proceeds by considering a variation 
Sip^it) = (p^it) — ip^it) between two paths (p^it), tp^it) (t G [ti, t 2 ]) with common ends: 



1 Note that Einstein's summation convention is extended to include also diagonal components of 
the inertia tensor. 
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</? M (£i,2) = ^(^1,2)- Together, they form a closed path in the space M^,. The above naive 
application of Hamilton's action principle in the system B employed analogous closed-path 
variations 5&{t) = &(t) - between two paths (t E [ti,fa]), with also 

common ends: $ l (ti,2) = 3>*(ti,2)- This, however, cannot be correct. Under a transformation 
(H), variations 5<^ M do not produce closed-path variations <5$ l in the space M$. The nonzero 
Burgers vector b l 7^ of the anholonomy causes a closure failure. This has to be accounted 
for in a correct derivation of the equations of motion. 

The anholonomy of the transformation (f|) requires distinguishing two types of variations 
of the paths in the space M$: closed path variations 5<& l (t) of the tpye described above, and 
anholonomic variations <5$ l (t), which are images of the 8^{t) variations in the space M^,. 
Only the latter will produce the correct equations of motion from an action principle in the 
system B. 

To calculate the variation of the classical action in the system B we derive the following 
simple formula for the anholonomic variation 5$* of the angular velocity Q l = 

5$* = is* + e^^W (15) 

(corresponding to formula (9) in Ref. [TT|] ). 

Indeed, the transformation @ and the definition (|7|) lead to the equation 

6& = 6(e\^n = e\5^ + d v e\^6<p v = j f (e\5^) + (d/ v - d v e\) <p v 8if? (16) 
which by (|TT|) becomes 

m = ^-5& + 2S ki i &5$ k (17) 
at 



leading to flI5|) via (|12|)- The variations 5^ are still performed in the system S. They can be 
mapped locally into holonomic closed-path variations 5$ J in the system B by <5$ l = e % . 
They have the fixed-endpoint property 5& l (ti^) = 0, reflecting the closed-path condition 
8^(^1,2) = in the system S. 



S 



Let us emphasize that our final nonholonomic variations (15) are completely intrinsic to 
the system B. Applying these to the action Ab, an integration by parts with the boundary 
conditions 5$*(ti,2) = leads to 

dA B = t dtli&3& = - [ t2 dtS&(t + OxLf (18) 
Jti Jti 

Since 5$* are now ordinary (holonomic) variations, we find the Euler equations (1). Thus 
the anholonomic action principle 

SA B = (19) 
produces the correct equations of motion in the body-fixed system B. 

4) Action principle for general motion in the body-fixed system B 

To extend the action principle to the general motion of a rigid body one usually chooses 
the point O in B to be the center of mass of the body and studies the movements of O and 
rotations around O. Then Ij are the body's principal moments of inertia. The configuration 
space is a 6-dimensional manifold M^ 6 ) = x SO (3) and the kinetic energy T = T trans +T rot 
consists of two terms: The first is the translational kinetic energy T trans = ^MV l V l in the 
system B the second is the rotational kinetic energy T rot = ^Iffifl 1 in the system B. 

In the system S, the position and orientation of the body are parametrized by the 
(holonomic) coordinates ^} 2 3 = {? M } M 12 6 j with x M being the mass center 
coordinates and ^ the Euler angles. 

The body-fixed basis vectors in B (see Fig. 1) are some function of the Euler angles 
£i(<p). Their components in the cartesian basis of the system S are ef((p). They form a 
3x3 orthogonal matrix s((p) = {s^), s^ 1 = (s 1 ^) = £ T ■ The components of the velocity of 
the center of mass in the system S are v M = if. The components of x M with respect to the 
basis Si are V 1 = e\x fl . 

The last relation permits us to introduce a new set of anholonomic coordinates X 1 de- 
scribing the center of mass motion in the system S as seen from the system B. By analogy 
with formula (0), we define the infinitesimal increments dX l as: 
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dX l = e\((p)dx^. (20) 

The set {X % , $*} i=1 23 = j^ 7 }^ 1 6 gi ves a complete set of an anholonomic coordinates for 
the configuration space R/ 3 - 1 x SO (3) of the moving body. 

Let us write down the anholonomic geometry of this space. Introducing 6x6 matrices, 
whose rows and columns are labeled by capital Latin and Greek letters, 



£ 



£'a), S = (Saz) = £ t £ 



(21) 



the geometry possesses the affine connection Tas A = £i A 9a£ i -e with the torsion tensor 
Sae A = T[ae] A 7^ and a vanishing Cartan curvature tensor -Ras.a" = 0, implying an 
nonvanishing Riemann curvature tensor Rat,, a" 7^ 0. This geometry is a simple extension 
of the 3-dimensional anholonomic geometry on SO(3) described in Section 2. It reflects the 
structure of the present configuration space Ry> x SO (3), i.e., the space of configurations of 
the rigid body comprising translations of the center of mass and rotations around it. 



As a direct consequence of the definition fl20[) , we have the relation V 1 = e l ^ = X 1 
which, together with (|7]), leads to the following simple form of the rigid body's Lagrangian 
in the system B: 



M ■ ■ ■ ■ 1 

L B = —X l X l + 
2 2 



(22) 



As before, a naive application of Hamiltons principle in the system B would produce a 
wrong equation MX 1 = for the motion of the center-of-mass. The correct anholonomic 
variation of the velocity 



d 



SX { = j f (SX l ) + e l jk (5&X k - &5X k ) . 



Indeed, the transformation (16!) and the definition 



d 



lead to 



dX l = 8{e\&) = j t (e\Sx^) + d x e\(x^ x - <f X ^). 



(23) 



(24) 



The variations 5x^ may be mapped locally into holonomic closed-path variations 5X l = 
e l fl Sx 11 in the system B with the property 8X l (t\fl) = reflecting the closed path con- 
dition 5x^(t lt2 ) = 0. Formula (23) follows using the equations x M = e^^X k , 5x fi = 
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e^ k 5X k {e^ k e k v = 6$), <p x = e*^', 5ip x = e X j5<5> j {e x j e j 1/ = 5 x )j, and the relation 
(dee T )ij = -e ijk d$ k . 

By applying the anholonomic variation (23) to the translational energy, we obtain 

5T trans = j^MX^X 1 ) - MiX* + :' lh A> J .X'- )<S.X : . (25) 

An integration by parts using of the fixed-ends conditions £X M (t 1)2 ) = leads to the following 
expression for the total anholonomic variation of the action Ab = J^ 2 Lsdt: 

SA B = f t2 dt(Ii&8$ + MX i dX i ) = (26) 

= - f t2 dt[6&{L + fl x L)* + 6X*(P + n x P)*] (27) 
Jti 

From the second term, we find the equations (2) for the translational motion. 
5) Conclusion 

By subjecting the action Ab in the body-fixed system to the new anholonomic variations 
with respect to translational and rotational degrees of freedom of the rigid body we have 
been able to derive both the correct Euler equations (|l|) and the equations (0) completely 
within the body-fixed system without reference to the stationary systems. 

The existence of such an action principle intrinsic to the body-fixed system may not be 
only of aesthetic value but may also have important practical consequences. 
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